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Abstract 

We develop a dimensional regularization approach to deal with the low-energy effects 
of the long-range Coulomb interaction in ID electron systems. The method allows us 
to avoid the infrared singularities arising from the long-range Coulomb interaction at 
D = \, providing at the same time insight about the fixed-points of the theory. We show 
that the effect of increasing the number of subbands at the Fermi level is opposite to 
that of approaching the bare Coulomb interaction in the limit — > 1. Then, we devise 
a double scaling limit, in which the large N effects are able to tame the singularities 
due to the long-range interaction. Thus, regular expressions can be obtained for all 
observables right at D = 1, bearing also a dependence on the doping level of the system. 
Our results imply a variation with in the value of the exponent for the tunneling 
density of states, which is in fair agreement with that observed in different transport 
experiments involving carbon nanotubes. As the doping level is increased in nanotubes 
of large radius and multi- walled nanotubes, we predict a significant reduction of order 
jV~^/^ in the critical exponent of the tunneling density of states. 
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1 Introduction 



In the last decade carbon-based devices have been largely studied as a possible new frontier 
of microelectronics. The limits of further miniaturization (predicted by Moore's law) have 
increased the research toward the development of molecular electronics and new efforts of 
scientists have been stimulated by the progress in carbon technology, in order to build a car- 
bon based microelectronics [1, 2]. Carbon nanotubes are the bricks of this new building and 
recent experiments have revealed that they are also excellent systems for the investigations 
of electronic transport in one dimension. It has been shown that the single-walled nanotubes 
can have semiconducting or metallic behavior, depending on the helical arrangement of the 
carbon rings around the tubule [3]. The metallic nanotubes have all a typical band struc- 
ture at low energies, in which two different subbands cross at opposite Fermi points (in the 
undoped system). 

It is well-known that the so-called Luttinger liquid behavior describes the regime with 
absence of electron quasiparticles that is characteristic of one-dimensional (ID) electron 
systems with dominant repulsive interactions [4, 5, 6]. In the case of the undoped metallic 
nanotubes, the fact of having four low-energy linear branches at the Fermi level introduces a 
number of different scattering channels, depending on the location of the electron modes near 
the Fermi points [7] . It has been shown, however, that processes which change the chirality of 
the modes, as well as processes with large momentum-transfer, are largely subdominant with 
respect to those between currents of like chirality [8, 9]. It has been concluded, therefore, that 
the metallic carbon nanotubes should fall into the Luttinger liquid universality class [9, 10]. 
Supporting that expectation, characteristic experimental signatures such as the power-law 
behavior of the tunneling density of states have been measured in ropes [11] as well as in 
individual single- walled nanotubes [12]. 

From the theoretical point of view, a relevant question is the determination of the effects 
of the long-range Coulomb interaction present in carbon nanotubes. In most part of the 
analyses carried out for such systems, the electron-electron interaction is taken actually as 
short-range, on the assumption that the finite length of the experimental samples has to 
impose in practice an infrared cutoff on the interaction. It is known, however, that the 
Coulomb interaction is not screened in one spatial dimension [13, 14], and that it remains 
long-range even after coupling a parallel array of ID systems — although the strength of the 
interaction can be consequently reduced in this case [15, 16]. 

The issue of considering the effects of the long-range Coulomb interaction is significant, 
as they have been shown to lead in general to unconventional electronic properties [17]. Such 
effects have been shown to be responsible for a strong attenuation of the quasiparticle weight 
in graphite [18]. The long-range Coulomb potential l^(|r|) ~ l/|r| yields the main electron 
interaction in carbon nanotubes as well as in the two-dimensional (2D) layers of graphite, 
which have a vanishing density of states at the Fermi level. An important confirmation of the 
marginal Fermi liquid behavior in such 2D layers comes from the experimental measurements 
in graphite of a quasiparticle decay rate linear in energy [19]. Owing to the singular nature 
of the Coulomb interaction, the imaginary part of the electron self-energy in the 2D layer 
behaves like (?^c<j at weak g coupling [20]. On the other hand, the effective coupling scales 
at low energy like g ~ l/log(cc;). This fact explains that the quasiparticle weight is not 
driven to zero by its logarithmic renormalization, being corrected instead by terms of order 
gHog{u:)^l/\og{u) [21]. 

The experiments on tunneling into the bulk of multi-walled nanotubes (MWNTs) have 
shown a power-law behavior of the tunneling density of states, with an exponent a ~ 0.3 
[22]. This is close to the values measured in single- walled nanotubes [11, 12], but it is 
much larger than expected in a wire with a large number N of subbands, where there 



should be a significant reduction in the strength of the electron-electron interaction. There 
have been attempts to confront this puzzle under the assumption of a disorder-enhanced 
renormalization of the interaction [23]. However, the exponents predicted in this treatment 
display a typical \/N reduction that make them to fall short with respect to the experimental 
estimate [22] . The enhanced suppression of tunnehng in the MWNTs supports therefore the 
view that the long-range character of the Coulomb interaction has to affect the system in a 
sizeable way. 

In a recent Letter [24] , we developed an analytic continuation in the number D of dimen- 
sions to accomplish the renormalization of the long-range Coulomb interaction at D — > 1. 
This limit turns out to be highly singular, since the marginal Fermi liquid behavior charac- 
teristic of the 2D graphene sheet is established as soon as one departs from D = 1. However, 
the attenuation of the electron quasiparticles is increasingly strong as D — 1, leading to an 
effective power-law behavior of observables like the density of states. In this way, we were 
able to predict a lower bound of the corresponding exponent that turned out to be very 
close to the value measured in experimental observations of the tunneling conductance for 
MWNTs [22]. 

In this paper, we focus on the singular behavior that arises at the limit D ^ 1 and 
introduce the efi^ect of the number of subbands that contribute to the low-energy properties 
of carbon nanotubes. This issue is relevant for the investigation of the nanotubes of large 
radius that are present in the MWNTs, which are usually doped and may have a large 
number of subbands crossing the Fermi level [25]. From a technical point of view, it becomes 
suitable the implementation of a large-A^ approach in the renormalization of the electron 
system. Increasing the number of subbands has the effect opposite to approaching the bare 
Coulomb interaction in the limit D ^ 1. We will see that it is possible to take a kind of 
double scaling limit, N ^ oo and — > 1, in which the large- A?" effects are able to tame 
the singularities related to the long-range interaction. Thus, regular expressions can be 
obtained for all observables right at D = 1, bearing also a dependence on the doping level 
of the nanotube system. 

The property of the double scaling limit, of rendering finite the effects of the renormalized 
long-range interaction, can be explained by the onset of some sort of dynamical screening, 
in the limit of an infinite number of subbands. Of course, in the real carbon nanotubes of 
a multi-walled sample, the number of subbands at the Fermi level may be large but finite, 
and the singularity of the Coulomb potential is not precisely reached due to the finite size 
of the system. However, the double scaling limit describes the universality class capturing 
the physics of the doped multi-walled samples, as any other ID approach to these systems 
is bound to bear the divergence of the long-range interaction in the low-energy limit. 

In this paper we deal with a renormalization group (RG) approach, in order to obtain 
the low-energy behavior of carbon nanotubes with long-range Coulomb interaction and a 
large number of subbands at the Fermi level. This approach is well-suited to the description 
of the scahng properties of observables like the quasiparticle weight or the density of states 
near the Fermi level. We will see that the renormalization introduces a factor of 1/ y/N in 
the critical exponent for those quantities. Our results can be compared with those from the 
experiments reported in Ref. [22], where measurements of the tunneling conductance have 
been carried out in doped MWNTs, with a number of subbands at the Fermi level iV pa 10 
(in the outer layer). For such values of N, and taking into account the long-range Coulomb 
interaction, the critical exponents that we obtain match well the values estimated from the 
experiments. 



2 Band structure of graphene and carbon nanotubes 



We describe the band structure of carbon nanotubes by the technique of projecting the band 
dispersion of a 2D graphite layer into the ID longitudinal dimension of the nanotube. The 
2D band dispersion of graphene can be found in Ref. [26]. It consists of an upper and a 
lower branch that only touch at the corners of the hexagonal Brillouin zone. Thus, when 
the system is at half-filling, the metallic properties derive from a pair of inequivalent Fermi 
points, around which there is conical dispersion for the modes of the graphene sheet. 

The 2D layers in graphite have a honeycomb structure with a simple hexagonal Br avals 
lattice and two carbon atoms in each primitive cell. We can start from the wavefunction 
associated to the electronic density, where we introduce two variational parameters, corre- 
sponding to two orbitals localized around the two carbon atoms in the primitive cell. Then, 
by introducing the wavefunctions in the energy functional for the lattice hamiltonian, we 
can solve the variational problem and obtain the following band dispersion, shown in Fig. 
1(a) [26]: 



E{k) = ±7^ 1 + 4cos2(^A;,) + 4cos(^A;,) cos(^y (2.1) 

Since the basis of the honeycomb lattice contains two atoms, there are two sublattices 
and two degenerate Bloch states at each Fermi point. If we choose the Bloch functions 
separately on each sublattice such that they vanish on the other, then we can expand the 
electron operator in terms of the Bloch waves 

*(,(x,y) ~ ^exp(-iQ;K • r)Fpaaix,y) (2.2) 

pa 

where a = ± labels the Fermi point, r = {x,y) lives on the sublattice p = ± under con- 
sideration and Fpaa{x,y) denote slowly varying operators. Thus, we can conclude that the 
low-energy excitations of the honeycomb lattice at half-filling are described by an effective 
theory of two 2D Dirac spinors [26] . 

Starting from the graphene band structure (2.1), after introducing periodic boundary 
conditions due to the cylindrical geometry of the tube, we obtain the energy bands of a 
carbon nanotube 




where N;, is the number of periods of the hexagonal lattice around the compact dimension 
(in the y direction) of the cylinder. So, at half-filling, metallic nanotubes have two Fermi 
points corresponding to large momenta ±Kf (see Fig. 1(b)). Now the low-energy expansion 
(2.2) transforms correspondingly and the electron operator is written as [27] 

-^Ax, y) = J2{27tR)''/^ exp(-mK • r) i/jp^Ax) (2.4) 

pa 

which introduces ID fermion operators ippaa{x) depending only on the longitudinal coordi- 
nate X. 

For later use, we note that Eq. (2.3) implies a simple dependence of the number of 
subbands that are found below the energy E = 

ns(e) — — arcsin(-\/l — £^) ~ —(^ + 77 + •••) 
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At low energies, we can substitute the function ns{s) with the first term of its Taylor expan- 
sion without a large disagreement. 



3 Dimensional regular izat ion near D = 



1 



In a previous paper [24], we have developed an analytic continuation in the number of dimen- 
sions in order to regularize the infrared singularity of the long-range Coulomb interaction at 
D — 1 [28]. Our aim was to find the effective interaction between the low-cncrgy modes of 
metallic nanotubes, which have linear branches crossing at the Fermi level. For this purpose, 
wc have dealt with the analytic continuation at general dimension D of the linear dispersion 
around each Fermi point. We start then with the hamiltonian 

/■A r 

B = "fE/^ dpipr-'/ (^V':»<r-pv„,(p) 

where the matrices are defined formally by {cTj, aj} = 2Sij. Here p(p) are density operators 
made of the electron modes ipaaip), and c{D)/\p\^^^ corresponds to the Fourier transform 
of the Coulomb potential in dimension D. Its usual logarithmic dependence on |p| at D = 1 
is obtained by taking the ID limit with c{D) = T{{D - I) /2) / {2^f-^ . 

A self-consistent solution of the low-energy effective theory has been found in [24] by 
determining the fixed-points of the RG transformations implemented by the reduction of 
the cutoff A. A phenomenological solution of the model was firstly obtained [29], carrying 
a dependence on the transverse scale needed to define the ID logarithmic potential, which 
led to scale-dependent critical exponents and prevented a proper scaling behavior of the 
model [29, 30]. The dimensional regularization approach of Ref. [24], which we follow here, 
overcomes the problem of introducing such an external parameter. 

In the present paper, we incorporate the effect of having a number of subbands crossing 
the Fermi level, which multiply consequently the number of electron fields and terms in the 
hamiltonian (3.1). Actually, as long as the long-range Coulomb potential is strongly peaked 
at p = 0, we will focus on the dominant interactions that take place between currents in 
which the electron modes remain in the same linear branch. 

Thus, the one-loop polarizability n(k, cj^t) is given by the sum of particle-hole contribu- 
tions within each branch. For a given subband, this leads in one spatial dimension to the 
well-known result [5] 

For a number A^ of different subbands, we take the analytic continuation to general dimension 
D [31] 

n(k, u,) = 2Nh{D) y^J^ll^l,.n),. (3-3) 

where h[D) = '^ ^^^^^^"^1)^(1^^-^'^^'^^ ■ In the case of carbon nanotubes we have N = 2, with 
the spin degeneracy already taken into account in (3.3). In the case of a thick nanotube 
with a large number of subbands, the doping modifies the number of subbands to be taken 
into account near the Fermi level. Moreover, it is also conceivable that, in the process of 
renormalization, A^ may be given effectively by a variable number of subbands depending on 
the energy cutoff, if the relevant energy scale of the problem is large enough, as we discuss 
in appendix. 

In the large- A^ limit we are interested in, the electron self-energy S(k, Uk) is obtained by 
dressing the bare Coulomb interaction with the polarization effects given by (3.3). We have 



then, at the dominant level, 



— ^ 

Go{k-p,UJk-l^p)jrrn^i — (3-4) 

where Go is the bare electron propagator. 

At general D, the self-energy (3.4) shows a logarithmic dependence on the cutoff at small 
frequency ujk and small momentum k. This is the signature of the renormalization of the 
electron field scale and the Fermi velocity. In the low-energy theory, with high-energy modes 
integrated out, the electron propagator becomes 

^ = -^^ Z~^{uJk-VFcr'k) 

where g = 2Nh{D)c{D)e' /v^, f{D) = ,n^,n,.u.\D/mD) ^ ^niD) = ^^3_^2v?+? - The 
quantity Z^/^ represents the scale of the bare electron field compared to that of the renor- 
malized electron field for which G is computed. 

The renormalizcd propagator G must be cutoff-independent, as it leads to observable 
quantities in the quantum theory. This condition is enforced by fixing the dependence of 
the effective parameters Z and Vi? on A as more states are integrated out from high-energy 
shells. We get the differential RG equations for the effective parameters g and Z 

A-^logZ(A) = V(-i)-^"+i 

^^^"^^ -MD)^-i{g) (3.6) 



n(3 - D) + 2 



oo 



The right-hand-side of Eq. (3.7) is a monotonous increasing function of g, for any dimen- 
sion between 1 and 2. Ai D — 2, the right-hand-sides of these equations can be summed up, 
and they coincide for = 2 with the functions obtained in the case of the undoped graphene 
sheet [21]. For D = 1 the function f3{g) vanishes, so that the ID model has formally a line 
of fixed-points as it happens in the case of short-range interaction. However, in the present 
model the variable g is sent to strong coupling in the limit D — > 1, and it remains to be 
checked the behavior of the RG flow in that regime. 

In order to approach the limit D — 1, we have to look for the asymptotic dependence 
on D of the functions appearing in the RG equations. We will see that this dependence 
appears as -D — 1 and D — 3 factors, revealing that these arc the two critical dimensions, 
corresponding to a marginal and a renormalizable theory, respectively. 



Starting with the function P{g), we need to carry out the sum at the right-hand-side of 
Eq. (3.7). This is given in terms of the hypergeometric special functions 

Since we need the /3 function near D = 1, we can approximate Boi^g) with the simple 
function Bi{g) 

To first order in D — 1, we get then 

As long as the coupling g depends implicitly on the dimension we also need to know 
its asymptotic behavior near D — 1. In practice, we can translate the analysis to the initial 
value of the couphng, g'o(-D) = 2Nh{D)c{D)e^ /vp- Carrying out an expansion near D — 1, 
we obtain 

where ({s) = J2kLi is the Riemann zeta function. So we can use in our calculations the 
simple expression 

Finally, it can be seen that the right-hand- side of Eq. (3.6) does not behave as a simple 
power law. This is due to the fact that the terms of the series vanish too slowly when n 
increases. The 7 function can be expressed in the form 

1(9) = ^^-^gTM (3.13) 
where the series T]:,{g) is available on tables. For D = 1 it has the simple expression 

T^{g) = ^^-^ (3.14) 

The scahng of the electron wave function near D = 1 is therefore given by 

7fe)-^y^(2-;rT^-^) (3.15) 

This coincides formally with the anomalous dimension that is found at D = 1 in the exact 
solution of the Luttinger model, what provides an independent check of our RG approach 
to the ID system. 

4 RG scaling and low-energy density of states 

Now we can rewrite the RG Eq. (3.7) near D — 1 as follows: 

ds ^ _/(£) 2(fl^ / _ 1\ 

dx 2N D ^ ^ \ y/TTg J ^ ' 



where x = — log(A). At D = 1 we have formally from Eq. (4.1) a line of fixed-points, which 
cover all values of e^/vp. In spite of the fact that g goes to strong coupling as D — * 1, 
we can see that the behavior in Eq. (3.12) does not suffice to cancel the D — 1 factor at 
the right-hand-side of Eq. (4.1). Therefore, we make sure in this way that the long-range 
Coulomb interaction is at a fixed-point of the RG a,t D — 1. 
We can also rewrite the RG equation (3.6) as follows: 

dlogjZ) ms-D 

- 2N 2 ^^^^^^ ^^-^^ 

Near D = 1 we find a crossover to a behavior with a sharp reduction of the quasiparticle 
weight in the low-energy limit A — > 0. All that is displayed in Fig. 5 in appendix, where 
we have represented the electron field scale square Z. If D is above ~ 1.2, we have a clear 
signature of quasiparticles in the nonzero value of Z at low energies, whereas for lower values 
of D the picture cannot be distinguished from that of a vanishing quasiparticle weight. 

Our target is to compare theoretical results with measurements of the tunneling density 
of states carried out in nanotubes of large radius [22] . The density of states n{e) computed 
at dimensions between 1 and 2 displays an effective power-law behaviour which is given by 

n{e) ^ 2NZ{e)\e\^-^ ^ 
\og{n{e)) ^ \og{2N)+\ogZ{e) + {D-l)\og{\e\) (4.3) 

The last term in Eq. (4.3) comes from the analytic continuation of the conical dispersion, 
which gives rise at D = 2 to the well-known behavior of the graphite layer, n{e) ~ \e\. 

In order to obtain the linear dependence of \og{n{e)) on log(|£|), we have to introduce 
the low-energy behavior of Z{e) in Eq. (4.3), by which 

log(n(£)) - {az - {D-l))x = aox (4.4) 

If we want to keep the linear behaviour of log(Z) in x — — log(|£|), we have to make an 
approximation. We can solve the RG equation (3.6) near a; = with a simple first order 
expansion, \og{Z) ~ ^{gQ)x, where Qq is the initial value of the coupling. Then we can analyze 
the dependence of az on the dimension D and at large N starting from the expression 

«z ^ ^^-^Tn{go)9o (4.5) 

4.1 Analytic continuation near D=l 

Following Ref . [24] , we can obtain a lower bound for the exponent of the density of states, 
for the purpose of comparison with experiments, by estimating the minimum of the absolute 
value of for dimensions ranging between D — 1 and D — 2. The evaluation can be 
carried out using Eq. (4.4), after plugging the expression of az in (4.5), with its strong 

dependence on D. 

In order to obtain an analytical expression for a^,, we can approximate T£,{g), starting 
from Ti{g), by using for 1 < D < 2 the simple approximation T£){g) ?a ^^^^Ti^g). Thus, we 
can give a simple formula for a^ rewriting Eq. (4.4) in the form 

.o^,oT.U^^^mR±^.^n-l) (4.6) 



We obtain a minimum value for |«/)| as a function of D when we introduce the expression 
of go{D) obtained in the previous section. From Fig. 2 we can see that the maximum value 



for an {o-m) corresponds to a dimension between 1 and 1.2. If the number of subbands is 
increased, then the value of I^mI decreases while the corresponding dimension approaches 1. 

Our analytic computation allows to understand the interplay between the number of 
subbands A'^ and the dimension D at which the minimum of the critical exponent |a£)| is 
obtained. We find that this is around \aB\ ~ 0.3 for not too large values of N, which is 
in agreement with the values measured experimentally. However, the fact that the minima 
are obtained above D = 1 stresses once again the singular character of the ID Coulomb 
interaction. The physical description can be clarified by noticing that the purely ID limit 
can be approached by increasing the number of subbands, as we propose below. 



4.2 Double scaling approach 

Here we exploit the dependence on the number of subbands as a way of regularizing 
the divergent effects of the Coulomb interaction, when D approaches the critical dimension 
D — 1. As it is easily seen inspecting the az function in Eq. (4.5), a large N can be 
combined with the divergence m D — 1 from the coupling constant in Eq. (3.12), so that 
both quantities regularize each other in a sort of double scaling limit. From the physical 
standpoint, this mechanism of regularization can be regarded as a kind of screening of the 
singular Coulomb interaction aX, D — 1. 

We use the above method of regularization near the critical dimension D = 1, in order 
to give an alternative evaluation of the critical exponent a. In the limit of large A'^, the 
dimension where the crossover effectively takes place, corresponding to the maximum value 
of ao-i approaches D = 1. In this regime, the last term in Eq. (4.6) vanishes and we can use 
the value of az as an estimate of a. 

In the hmit of large A^ and small values of — 1, we have Ti{g) — > — y^, so that Eq. 

(4.5) becomes 

/(L')(3 -£))(£> + 1) , 

Taking into account the form of in Eq. (3.12), we see that az vanishes at large A^ as 
~ 1/ \/iV, and it diverges at D — >■ 1 as ~ 1/ \/ D — 1. 

Now we can use the regularization from the number of subbands, in order to access a 
physical regime that would be otherwise hidden in the usual treatment of the ID Coulomb 
interaction. In other words, we take simultaneously the limits A^ ^ oo and D ^ 1 by 
keeping a finite renormalized value A^^en for the number of subbands, A^ren = A^(-D — 1). We 
use then A^^en as a way of parameterizing the dependence of the critical exponent a on the 
doping level, which sets the number of subbands at the Fermi level in the carbon nanotubes. 
From Eq. (4.7) we obtain a simple expression for a as a function of the renormalized number 
of subbands 

The agreement between the two estimates of the critical exponent a by the double scaling 
limit and by analytic continuation near D = 1 is shown in Fig. 3(a). The dependence of az 
on the number of subbands in the second approach can be actually borne out by numerical 
results. These are in good agreement with our analytical predictions, as we show in Fig. 
4. The dependence of az on the number of subbands is rather similar to l/\/iV, especially 
for dimensions near 1. Thus, from this overall picture we can predict that, by changing 
the doping level in carbon nanotubes of large radius, a sizeable variation of the anomalous 
electron dimension should be observed in the range A^ ~ 2 — 100, with a consequent change 
in the power- law behavior of the tunneling density of states given by Eq. (4.8). 



5 Discussion 



In this paper wc have described the appropriate framework to deal with the low-energy 
effects of the long-range Coulomb interaction in ID electron systems. Our interest has been 
focused in the scaling behavior of quantities like the quasiparticle weight or the low-energy 
density of states, which can be compared directly with the results of transport experiments. 
For this purpose, we have developed a RG approach in a dimcnsionally regularized theory 
devised to interpolate between the 2D conical band dispersion of graphene and the ID band 
dispersion of the carbon nanotubes. The analytic continuation in the number of dimensions 
has allowed to avoid the infrared singularities that the long-range Coulomb interaction poses 
at D = 1, providing insight at the same time about the fixed-points and universality classes 
of the theory in the limit D ^ 1. 

Thus, we have seen that the ID electron system with long-range Coulomb interaction is 
at a fixed-point of the RG flow, for arbitrary values of the interaction strength. However, 
the scaling behavior of observables like the quasiparticle weight is dictated by divergent 
critical exponents in the limit D ^ 1. In order to make contact with the results of transport 
experiments on MWNTs, we have introduced the relevant dependence on the number 
of subbands at the Fermi level. By increasing the value of N, it has been possible to 
approach the critical dimension D — 1 while producing finite, sensible estimates of the critical 
exponents in the carbon nanotubes. The double scaling limit, N ^ oo and gives 
rise to a phase of the system that captures the physical properties measured in nanotubes of 
large radius, which use to be significantly doped. In this way, we have been able to describe 
a new universality class that arises in the large- A'^ limit, in which the effect of the large 
number of subbands is to regularize the singular ID Coulomb interaction. The new phase 
arises from a sort of ID screening of the long-range interaction, giving rise to a modified 
Luttinger liquid picture in which the critical exponents depend on the doping level as well 
as on the strength of the interaction. 

The framework that we have introduced is appropriate for the description of carbon 
nanotubes of large radius, in which the shift of the Fermi level leads to a sizeable variation 
of the number of subbands contributing to the low-energy electronic properties. Such a 
situation is realized in the MWNTs, which are concentric arrangements of carbon nanotubes 
with a diameter of ~ 10 nm and typically ~ 5 inner cylinders. The transport measurements 
carried out in the MWNTs refiect usually the electronic properties of the outer layer, to 
which the electrodes are attached. In these systems the number of subbands at the Fermi 
level can be large, i.e. of the order of ~ 10 — 20 taking into account the electron spin, because 
of the electrochemical doping. 

One of the most significant observations made in the MWNTs has been the power-law 
behavior of the tunneling conductance as a function of the temperature or the bias voltage. 
In the experimental measurements, the maximum energy scale reached in those variables has 
been always below 100 meV. This is of the order of the typical energy spacing between the 
different subbands, so that the transport measurements refer to the contribution of a fixed 
number of subbands at low energies. In the case that the experiments were carried out over 
a larger energy range, the physical description should be modified to allow for a variable 
number of subbands depending on the cutoff energy, producing the results that are shown 
in appendix. 

The measurements carried out in the MWNTs have displayed a power-law behavior of 
the tunneling conductance, that gives a measure of the low-energy density of states, with 
exponents ranging from 0.24 to 0.37 [22]. These values are, on the average, below those 
measured in the single- walled nanotubes, which are typically about 0.35 [12]. We see 
that our results can account satisfactorily for this slight reduction in the critical exponent 



with the change of the nanotube thickness. We observe from Fig. 3(b) that, for a number 
of subbands N = 2, the vahie of |a| corresponding to an imdoped single- walled nanotube 
is slightly below 0.3. This means that our choice for the interaction strength (with some 
reduction from the dielectric constant) has been appropriate. The important point is that, by 
varying the number N of subbands from 2 to 10, we obtain a shift in the value of the critical 
exponent which is in fair agreement, on the average, with that observed from the different 
transport experiments. This lends strong support to our picture of the renormalization of 
the long-range Coulomb interaction. 

The main prediction that comes from our study is that there should be a significant 
reduction in the critical exponent of the tunneling density of states as the doping level is 
increased in nanotubcs of large radius. Wc have seen that the presence of a large number of 
subbands at the Fermi level implies a reduction of \a\, by a factor of the order of N'^^"^. It 
would be of considerable importance to test such a dependence in experiments carried out 
using samples with different amount of doping. That kind of effect should arise also as a 
consequence of considering wider compact nanotube structures, stressing in another way the 
significance of the geometry in the electronic properties of the nanotubes. 
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A Number of subbands depending on cut-ofF energy 



Next, we consider a number of subbands A''^ depending on the cut-off. So, we introduce tiie 
energy dependent function ns{e) which substitutes the usual Us- We define a unit energy 
as Sp^ so that we have Xp — — log(£p) = 0. The corresponding number of subbands reads 
N,{ep) ^No + Np [32] 

Ns{e) = No + Np-^ 

Ep 

Ns{x) ^ No + NpC-'' ^ = 1 + npe''' 

Here Up represents the number of subbands, in a unit of energy normahzed with respect to 

When Hg depends on the energy we obtain, by increasing the number of subbands, that 
az vanishes as However a stronger contribution causes a growth of a (although 

its absolute value stays less than 1). The relevant correction in the eq. (4.4) is due to the 
direct contribution of the number of subbands depending on the energy. This is the ccat term 
which can strongly modify the ao value. Near a; = 0, we can obtain this term as a simple 
function of rip, which represents the number of subbands, in a unit of energy normalized 
with respect to A^o 



This term does not depend on the dimension, so it does not affect the search of the maximum 
in the q;/5, rather it gives a simple shift in the global a value, which we have to compare 
with the experimental value of the parameter. Larger values of Up give larger contributions 
(in absolute value) from this term which ranges between and 1 (once again, in absolute 
value) . 

In the case of a linear dependence of the number of subbands on the energy, a large 
correction to the measurable a exponent can arise from the contribution (oat ~ — 1), 
while the effective parameter Z yields strong deviations from the behaviour at n = 4, even 
when the crossover dimension approaches the fixed point D = 1. 

Some interesting effects regarding the numerical calculation of Z arise, when the number 
of subbands depends on the cutoff energy {rip is the density of the number of subbands in 
the energy unit). In this case a strong deviation from the usual behaviour appears. The 
renormalization of the electron field scale Z . for a fixed number of subbands (4 subbands, 
i.e. iV = 2), shows a drastic suppression of the electron field scale, which takes place over a 
variation of only two orders of magnitude in the energy scale. When the number of subbands 
increases, we observe that the scale of energies, and also the shape of the curve, are quite 
different. The energy scale goes towards lower and lower energies, while a large plateau 
appears in the curve Z{—\og[E)). A simple mechanism can explain this behaviour. When 
the cutoff increases, many subbands have to be taken into account and the system, at any 
dimension greater than 1, is more similar to a 2-dimensional graphene sheet than to a 1- 
dimensional wire. Only when the energy goes below a small value (i.e. x > xt) the system 
restores its Luttinger liquid behaviour and the quasiparticle amplitude vanishes, following 
the usual behaviour. This is what we can observe on the left of Fig. 5 . All that is clear 
also on the right, where we show the dependence of az ~ ''^^^dx^^ Here the asymptotic 

behaviour reflects the small Up results, whereas the values for 'high energies' are quite similar 
to those of the usual 2D graphene. 

All the phenomenology of a tube with a very large number of subbands below the cutoff 
energy is more similar to that of 2D graphene than to a ID Luttinger liquid. This effect 
could be seen either as a smaller crossover dimension, or as the persistence of quasiparticle 



poles at relatively high energy. We can obtain the same physics if we imagine the radius of 
the tube going to infinity: the geometrical limit of the tube is the graphene layer, while the 
energy spacing between two subbands vanishes, so that an infinite number of subbands have 
to be taken into account, also for very small energies. 
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Figure Captions 



Fig. 1: (a) Energy levels in a graphite sheet: the cusps appear at the six corners of the first 
Brillouin zone, (b) Low-energy band structure of metallic carbon nanotubes with different 
radius. The energy is measured in units of the hopping parameter and the momentum in 
units of the inverse lattice spacing. As shown in (c), the radius is connected to the value of 
Nb in a simple way: R = ^/Sa / {2% sm{7r / Nh)) ^ A^6\/3a/(27r), where a denotes the honey- 
comb lattice constant (a/\/3 = 1.42A). 

Fig. 2: In (a) we show a as a function of the dimension (Eq. (4.6)), and we choose 
the maximum value of this function between D = 1 and D = 2 as an estimate of the critical 
exponent. We find a value for a in the usual N = 2 undoped nanotube that reproduces the 
anomalous exponent measured experimentally {a ~ —0.3), corresponding to a dimension for 
the crossover between 1.1 and 1.2. The value of a is compared to the term az{D) (dashed 
line) dominant near the divergence — > 1 . As shown in (b) , when the number of subbands 
increases, the absolute value of the maximum decreases, while the corresponding dimension 
approaches 1 (each curve corresponds to a different value of N) . 

Fig. 3: In this paper we propose also a different way to estimate the critical exponent, 
as the value of az in the double limit N ^ oo and D — > 1. A comparison between the re- 
sults of the double scaling approach (dashed line) and those from the analytic continuation 
in D (full line) is shown in (a). The increase of a for a moderate number of subbands is 
shown in (b). We have taken in all cases a value of the coupling /{ti'^vf) ~ 0.5, which is 
appropriate for MWNTs, owing to the reduction in this value due to the interaction with 
the inner metallic cylinders [15]. We can conclude that, when N is very large, there is a 
good agreement between the results obtained with the two different approaches, i.e. double 
scaling and analytic continuation. The dashed portion of the curve in (b) agrees, in the 
one-mode limit (i.e. for N = 0.5), with the results of Ref. [24]. 

Fig. 4: On the left, we show the linear dependence of az on the inverse square root of 
N, for different values of D. On the right, wc observe a rather smooth behavior of the slope 
of the curves for each dimension in the first figure. 

Fig. 5: On the left, we plot Z{—\og{E)) for a different number of subbands, depend- 
ing on the energy. On the right, we show the dependence of 
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